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ABSTRACT 

We present detailed NLTE spectral synthesis models of the Type II supernova 
2005cs, which occurred in M51 and for which the explosion time is well deter- 
mined. We show that previous estimates for the reddening were significantly 
too high and briefly discuss how this will effect the inferred progenitor mass. 
We also show that standard CNO-burning enhanced abundances require far too 
large an oxygen depletion, although there is evidence for a single optical N II 
line and the sodium abundance shows clear evidence for enhancement over solar 
both as expected from CNO processing. Finally we calculate a distance using 
the SEAM method. Given the broad range of distances to M51 in the literature, 
the determination of a distance using Cepheid variables would be quite valuable. 

Subject headings: stars: abundances, atmospheres — supernovae: SN 2005cs 



1. Introduction 



Supernova 1987A (SN 1987A) confirmed our basic theoretical understanding of a Type II 
supernova (SN II) as the core collapse of a massive star, which leaves behind a compact object 
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(neutron star or black hole) and expels the outer mantle and envelope into the interstellar 
medium (see lArnett et al.l Il989l . and references therein). While the explosion mechanism 
itself remains a subject of active research, the propagation of the shock wave through the 
mantle and envelope is reasonably well underst ood. Theoretica l models for the light curve 
do a good job of reproducing the observations (jBlinnikovl If 9991 ; iBlinnikov et al.l |2000| ) and 
detailed radiation transport calc ulations verify that these m odels reproduce the observed 
spectra from the UV to the IR ([Mitchell et al.l l200ll . 120021 ) . Unlike Type la supernovae 
(SNe la) where the actual compositions as a function o f velocity are an important subject 
of current research (IFisher et al.lll999l ; IStehle et al.ll2005l ). the compositions of the envelopes 
of red or blue supergiant stars is primarily hydrogen and helium and thus we show that the 
primordial compositions can be determined by the detailed spectral modeling of observed 
SN II spectra. Of particular interest is the effect of CNO processing on the composition of 
the ejected material. 

SNe II have a very large spread in their intrinsic brightness, from the very dim SN 1987A, 
to the exceedingly bright SNe 1979C and 1997cy. The observed spread in intrinsic luminosity 
is greater than a factor of 500. This is not surprising given the fact that the progenitors span 
a wide range of initial stellar masses, possible binary membership, and prior star formation 
histories. Clearly, SNe II do not meet the astronomers requirement of being a standard 
candle; however, since our models determine the stellar compositions and total reddening 
this provides evidence that their atmospheres can be well understood. These results in- 
crease the attractiveness of SNe II as cosmological probes and in the future SNe II will 
become complementary with SNe la as dis tance indicators. Both the SEAM (Spe ctral- 



fitting Expanding Atmosphere Method) (see lBaron et al.lll995bl . Il996l ; iMitchel 



et al J 12002. 
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tances to SNe II (and other supernova types for SEAM) depend on the ability to model the 
spectral energy distribution (SED) of the supernova atmosphere accurately. Clearly, quan- 
tities like the composition and total reddening play a role in the output SED and thus the 
work presented here helps to place both these methods on firmer footing. 

While the average SN II has a luminosity several times less than a typical SN la, cur- 
rent ground-based searches and proposed future space-based searches for supernovae will 
easily detect these objects at cosmologically interesting distances. Using an empirical re- 
lation between brig htness and expansion velocity on the plateau portion of the light curve 
Nugent et al.l (120061 ) were able to construct a Hubble diagram using 5 SNe II out to a redshift 
z = 0.3. 
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The explosion mec hanism is thought to result in a non-spherical ly symmetric shock- 
wave (see for example iLeBlanc &: Wilson! Il970l ; iKhokhlov et al.l Il999l ) and t his has been 



confi r med by the detection of significant polari zation in the spectra of SNe II (I Wang et al. 



19961 ; iLeonard et al.l |2000| ; iLeonard et al.l 120061 ) . Nevertheless, the large expans ion of the 



eject a before the supernova becomes visible leads to significant "sphericalization" (j Chevalier 
19841 ). and asphericity effects may only be important at late times, or in SNe II that are 
strong circumstellar interacters (SNe Iln). These supernovae are observationally clearly 
distinguishable from "normal" Type II supernovae. 



Calculations (IMiralda-Escude fc Reeall997l ; lDahlen &: Franssonlll999l ; lYungelson fc Livio 



20001 ) indicate that SNe II may in fact be the first stellar objects visible in the universe and 
hence they serve as important probes for the star formation rate, the rate of chemical evolu- 
tion by measuring their primordial abundances, and directly for the cosmological parameters. 
We show that the primordial abundances of the progenitor star can be determined by detailed 
synthetic spectral modeling. In this work we do not seek to determine exact nucleosynthetic 
yields, but rather show that this is feasible given a large grid of models that will become 
available as more supernovae are analyzed and as computational resources become more 
plentiful. 



2. Calculations 



We have chosen to model SN 2005cs in detail because there are several high signal to 
noise spectra and UBVRI photometry very near the time of explosion as w ell as at later 



times . SN 2005cs was discovered in M51 (NGC 5194) on 2005 Jun e 28.905 (IKloehr et al. 
2005h . the earliest detection was that of M. Fiedler on 200 5 June 27.91 JPastorello et al-lbood ) 



and nothing was visible on June 20.6 (IKloehr et al.ll2005f). Moreover no other detections were 
found on June 26 by other amateur observers (IPastorello et al.ll20061) . thus the explosion d ate 
is know n to within abou t 1 day . Using the rather long distance of iFeldmeier et al.l (119971 . see 
below). IPastorello et al.l (120061 ) classified SN 2005cs as moderately under-luminous compared 
to typical SNe IIP and in the clas s of SN 1997D and SN 1999br. In addition a progenitor 
has been identified for SN 2005cs (IMaund et al.ll2005l ; iLi et al.l 120061 ) with a mass of about 
10 M . 

The calculation s were performed using the multi-purpose stellar atmospheres program 



PHOEN IX versio n 14 (lHauschildt fc Baronlll999l : Baron fc Hauschildtlll998l : lHauschildt et al. 



1997al lbl. Il996h . PHOENIX solves the radiative transfer equation along characteristic rays 
in spherical symmetry including all special relativistic effects. The non-LTE (NLTE) rate 
equations for many ionization states are solved including the effects of ionization due to 
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non-thermal electrons from the 7-rays produced by the radiative decay of 56 Ni, which is 
produced in the supernova explosion. The atoms and ions calculated in NLTE are: H I, 
He I II, C I-III, N I-III, O I-III, Ne I, Na I-II, Mg I-III, Si I III, S I III, Ca II, Ti II, Fe I III, 
Ni I-III, and Co II. These are all the elements whose features make important contributions 
to the observed spectral features in SNe II. 

Each model atom includes primary NLTE transitions, which are used to calculate the 
level populations and opacity, and weaker secondary LTE transitions which are included in 
the opacity and implicitly affect the rate eq uations via their effect on the solution to the trans- 



port equation (lHauschildt fe Baron! Il999l ). In addition to the NLTE transitions, all other 
LTE line opacities for atomic species not treated in NLTE are treated with the equivalent 
two-level atom source function, using a thermalization parameter, a = 0.05. The atmo- 
spheres are iterated to energy balance in the co-moving frame; while we neglect the explicit 
effects of time dependence in the radiation transport equation, we do implicitly include these 
effects, via explicitly including the rate of gamma-ray deposition in the generalized equation 
of radiative equilibrium and in the rate equations for the NLTE populations. 

The models are parameterized by the time since explosion and the velocity (vq) where 
the continuum optical depth in extinction at 5000 A (r st d) is unity, which along with the 
density profile determines the radii. This follows since the explosion becomes homologous 
(v oc r) quickly after the shock wave traverses the entire star. The density profile is taken 
to be a power-law in radius: 

p oc r~ n 

where n typically is in the range 6 — 12. Since we are only modeling the outer atmosphere of 
the su pernova, this simple pa rameterization agrees well with detailed simulations of the light 



20001 ) for the relatively small regions of the ejecta that our models 



curve (IBlinnikov et al. 
probe. 

Further fitting parameters are the model temperature T mo d e i, which is a convenient way 
of parameterizing the total luminosity in the observer's frame. We treat the 7-ray deposition 
in a simple parameterized way, which allows us to include the effects of nickel mixing which is 
seen in nearly all SNe II. Detailed fitting of a time series of the observed spectra determines 
all the parameters. 



3. Reddening 



Determining the extinction to SNe II is difficult, since they are such a heterogeneous 
class, it is difficult to find an intrinsic feature in the spectrum or light curve that can be 
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used to find the parent galaxy extinction. iBaron et al.l (120031 . 120001 ) found that the Ca 11 
H+K lines can be used as a temperature indicator in modeling very early observed spectra. 
This suggests that detailed modeling of SNe IIP spectra can be used to determine the total 
reddening to the supernova. F or SN 2005cs the reddening has been estimated in a variety 
of ways. iMaund et al.l (120051 ) used the relationship between the equivalent width of the 
Na I D interstellar absorption line to obtain a color excess of E(B — V) = 0.16, as well 
as the color magnitude diagram of red supergiants within 2 arcsec of SN 2 005cs to obtain 
E(B - V) = 0.12, and their final adopted E(B - V) = 0.14. Ili et al.l tood ) noted the large 
scatter in the relationships for equivalent width of the Na I D line, obtaining a range of 
E(B — V) — 0.05 — 0.20. Assuming that the color evolution of SN 20 05cs is similar to t hat o f 
SN 1999em, they found E(B - V) = 0.12. Also using the Na I D line IPastorello et all (hood ) 
found E(B — V) = 0.06, but noting the uncertainty and comparing with the work of other 
authors the y adopted E(B — V) = 0.11. We began our work by adopting the reddening 
estimate of IPastorello et al.l (120061 ) . since our spectra were obta ined from these authors . 
Figure [1] (top panel) shows our best fit using solar abundances (IGrevesse fc Sauvall 11 998) 



where the observed spectrum has been dereddened using the reddening law of ICardelli et al. 
(119891 ) and Ry = 3.1. Figured] (bottom panel) shows the same for the case of CNO enhanced 
abundances ( Dessart Sz Hillier 2005 : Prantzos et al.l[l986l ). From the figure it is evident that 
the region around H/3 is very poorly fit, there is a strong feature just to the blue of H/3 and 
H/3 itself is far too weak. We attempted to alter the model in a number of ways, changing the 
density profile, velocity at the photosphere, and gamma-ray deposition in order to strengthen 
H/3, however we were unable to find any set of parameters that would provide a good fit to 
H/3 (and the rest of the observed spectrum) with this choice of reddening. These models have 
T mo dei=18000 K, vo— 6000 km s _1 , and n = 8. After careful analysis the strong feature too 
the blue of H/3 is due to C III and N III lines that become strong at these high temperatures. 
Figure [2] shows that if T model is reduced to 1200 K and the color ex cess is reduced to the 
galactic foreground value of E(B — V) — 0.035 (ISchlegel et al.lll998l ) the fit is significantly 
improved. Our yalue of E(B — V) = 0.035 is in agreemen t within the errors of the lower 
values found by iLi et al.l (120061 ) and IPastorello et al.l (120061 ) . Clearly, the Na D interstellar 
absorption indicates that there was some dust in the host galaxy, but we have chosen the value 
that gave the best fit for this epoch which was no dust extinction in the host galaxy. Given 
our modeling uncertainties, E(B — V) — 0.035 — 0.05, but we will adopt the value of 0.035 
for the rest of this work. Note that this lowe r value of the extin ction will som e what lower 
the inferred mass of the progenitor found by IMaund et al.l (120051 ) and ILi et al.l (120061 ) , but 
other uncertainties such as distance and progenitor metallicity also play important roles in 
the uncertainty of the progenitor mass. Comparing with the published results, the difference 
due to our distance and reddening estimates on the value of the pr ogenitor mass obtained 
is quite small. Using the value of Ry — Ri = 1.28 and Ri = 1.89 (ITammann et al.l 120031 ) 
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the change in / is ~ 0.15 mag and ( V — I) ~ 0.13 mag. Our derived distance is consist ant 
with t hat adopted by lLi et al.l (120061 ) and only somewhat smaller than that of iMaund et al. 
J2005h . 



Clearly the bluest part of the continuum is better fit with the T mo d e i=18000 K models 
than with the T mo d e i= 12000 K models. We did not attempt to fine tune our results to 
perfectly fit the bluest part of the continuum since the flux calibration at the spectral edges 
is difficult and it represents our uncertainty in T mode i and E(B — V). None of the above 
models have He I A5876 strong enough. It is well known that Rayleigh- Taylor instabilities 
lead to mixing between the hydrogen and helium shells thus our helium abundance is almost 
certainly too low, but we will not explore helium mixing further in this work. 
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Fig. 1. — . The observed spectrum on Jul 02, 2005 is compared to a model spectrum using 
solar abundances (top panel) and CNO enhanced abundances (bottom panel). The observed 
spectrum has been dereddened using E(B — V) — 0.11. In this and subsequent plots the 
model and observations have been normalized at 8000 A. 




Fig. 2. — . The observed spectrum on Jul 02, 2005 is compared to a model spectrum using 
solar abundances (top panel) and CNO enhanced abundances (bottom panel). The observed 
spectrum has been dereddened using E(B — V) = 0.035. 
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Abundances 



Typically, the approach to obtaining abundances in differentially expanding flows has 
bee n through line identificati ons. This method has been successful using the SYNOW code 

and referenc es therein) as well as the work of Mazzali and col- 



Branch et al. 2005 



2006 



Stehle et al.ll2005l ). Nevertheless, line identifications do not provide 



(see 

laborators (for example 
direct information on abundances, which are what are input into stellar evolution calculations 
and output from hydrodynamical calculations of supernova explosions and nucleosynthesis. 
Line identifications are subject to error in that there may be another candidate line that is 
not considered in the analysis or there may be two equally valid possible identifications, the 
classic being He I A5876 and Na I D, which have very similar rest wavelengths and are both 
expected in supernovae (and have both been identified in supernovae). 

Here we focus on nitrogen and sodium. Massive stars which are the progenitors of SNe II 
are expected to undergo CNO processing at the base of the hydrogen envelope followed by 
mixing due to dredge up and meridional circulation. This would lead to enhanced nitrogen 
and sodium and dep leted carbon and oxygen . For our CNO processed models we take the 
abundances used by iDessart fc Hillierl (120051 ). However, significant mixing of the hydrogen 
and helium envelope is expected to occur during the explosion due to Rayleigh- Taylor insta- 
bilities and mass loss will occur during the pre-supernova evolution. Studies of the evolution 
of massive stars show that for an 11 M model the mass fractions of CNO are very close 
to solar, with significantly less depletion of carbon and oxygen than is obtained from CNO 
processing (A. Chieffi & M. Tavani, 2006, in preparation). A detailed study of SN 2005cs 
spectra using the results of these hydro models will be a sensitive test of the predictions of 
stellar evolution theory. 



4.1. Line IDs 

In detailed line-blanketed models such as the ones presented here line identifications 
are difficult since nearly every feature in the model spectrum is a blend of many individual 
weak and strong lines. Nevertheless, it is useful to attempt to understand just what species 
are contributing to the variations in the spectra. In order to do this we produce "single 
element spectra" where we calculate the synthetic spectrum (holding the temperature and 
density structure fixed) but turning off all line opacity except for that of a given species. 
Figure [5] shows the single element spectrum for N II for our T mo d e i = 12000 K models with 
solar abundances (top panel) and CNO enhanced abundances (bottom panel). Clearly the 
N II lines are present in CNO enhanced models and absent in the solar models, but their 
effect on the total spectrum is unclear. 
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N II lines were f i rst id entified in SN II in SN 1990E JSchmidt et al.l Il993f ). Using 
SYNOW, iBaron et all (120001 ) found evidence for N II in SN 1999em, however more detailed 
modeling with PHOENIX indicated that the lines were in fact due to high velocity Balmer 
and He I lines. Promin ent N II lines in the optical are N II A4623, A5029, and A5679. 



Dessart fc Hillieri (120051 1 found strong evidence for N II in SN 1999em. Using SYNOW 



Elmha mdi et al. (in preparation) found evidence for N II in SN 2005cs, as did lPastorello et al. 



(120061 ) using the code developed by Mazzali and colla borators. Figure [2] compares solar abun- 



dances to a model with enhanced CNO abundances ( Dessart fc H illier 2005: iPrantzos et al. 



19861 ). Figure [2] shows that the CNO enhanced abundances model does somewhat better 



fitting the emission peak of H/3 (see below), however the feature to the blue of H/3, clearly 
well-fit in the solar abundance model is completely absent in the model with enhanced CNO 
abundances. Figure [3] shows the region where the optical N II lines are prominent and in 
particular, the N II A5679 is not quite in the same place as the observed feature and the two 
bluer lines have almost no effect. 
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Fig. 3. — The observed spectrum on Jul 02, 2005 is compared to a model spectrum using 
CNO enhanced abundances and one using solar abundances. The region of N II lines is 
shown for clarity. The vertical dashed lines show the rest wavelength of the N II optical 
lines, while the solid vertical lines show those same lines blueshifted by vq = 6000 km s _1 . 
The observed spectrum has been dereddened using E(B — V) = 0.035. 
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Fig. 4. — The observed spectrum on Jul 02, 2005 is compared to a model spectrum using 
CNO enhanced abundances and two where the oxygen abundance is enhanced (over the CNO 
value by a factor of 10 and nearly 100) in an attempt to improve the H/3 emission feature 
and still retain the desired O II feature just blueward of H/?. The plot shows how complex 
line IDs are. The alteration of the oxygen abundance not only affects the features, but the 
structure of the entire atmosphere itself (since oxygen is such a good electron donor). Neither 
modification is satisfactory and clearly more work needs to be done to obtain the "right" 
CNO abundances. The vertical dashed lines show the rest wavelength of the N II optical 
lines, while the solid vertical lines show those same lines blueshifted by vq = 6000 km s . 
The observed spectrum has been dereddened using E(B — V) = 0.035. 
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We attempted to further improve the fit to the H/3 emission feature by using the same 
abundances for C, N, and Na as our standard CNO processed abundances, but enhancing 
oxygen by a factor of 10 and nearly 100 (i.e. using the solar value in the latter case). Figure H] 
shows how complex line IDs are. The alteration of the oxygen abundance not only affects 
the features, but the structure of the entire atmosphere itself (since oxygen is such a good 
electron donor). Neither modification is satisfactory and clearly more work needs to be done 
to obtain the "right" CNO abundances. Interestingly, the model with enhanced nitrogen 
and solar oxygen clearly shows the N II A4623 line. 

In an attempt to identify the better fit of the feature just to the blue of H/3 we examined 
the single element spectrum of C II (Fig. [6]) since we considered that the C II A4619 line 
might be playing a role. Again this feature is clearly present in the solar model and absent 
in the CNO enhanced model. 

Finally, we examined the lines produced by O II. Figure [7] clearly shows that O II 
lines play an important role in producing the observed feature just blueward of H/?. Most 
likely it is the lines O II A4651.5 and A4698 which are producing the observed feature. On 
the other hand it is also clear that O II A4915 and A4943 are producing the deleterious 
feature just to the red of H/3. Since we want somewhat higher oxygen abu ndances, it was 



natura l to examine a model using the recently revised solar abundances of lAsplund et al. 



( 120051 ). Compared with our standard solar model abundances there was not a significant 



improvement (nor a degradation) in the fit to favor one over the other. 

Figure [8] examines the possibility that both C II and O II contribute to the feature just 
blueward of H/?. It is obvious that only O II lines are important for the structure of this 
feature. 

Thus, it is clear that the strong depletion of oxygen expected from CNO processing is not 
evident, we can not rule out that there is some enhanced nitrogen in the observed spectra, 
but the N II lines don't seem to form in the right place. However since we are studying 
simple, parameterized, homogeneous models this could be an artifact of our parameterization. 
Nevertheless, the absorption trough of the feature that we would like to attribute to N II 
A5679 is too fast in our models, whereas one would expect the CNO processed material with 
higher nitrogen abundance at the outermost part of the envelope and thus to form at even 
higher velocity due to homologous expansion. Clumping could of course change this simple 
one-dimensional picture. 
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Fig. 5. — The observed spectrum on Jul 02, 2005 is compared to a model spectrum where only 
N II line opacity has been included. The top panel is for solar abundances and the bottom 
panel using CNO enhanced abundances. The vertical dashed lines show the rest wavelength 
of the N II optical lines, while the solid vertical lines show those same lines blueshifted by 
Vq = 6000 km s™ 1 . The observed spectrum has been dereddened using E(B — V) = 0.035. 
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Fig. 6. — The observed spectrum on Jul 02, 2005 is compared to a model spectrum where only 
C II line opacity has been included. The top panel is for solar abundances and the bottom 
panel using CNO enhanced abundances. The observed spectrum has been dereddened using 
E(B - V) = 0.035. 
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Fig. 7. — The observed spectrum on Jul 02, 2005 is compared to a model spectrum where only 
O II line opacity has been included. The top panel is for solar abundances and the bottom 
panel using CNO enhanced abundances. The observed spectrum has been dereddened using 
E(B - V) = 0.035. 
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Fig. 8. — The observed spectrum on Jul 02, 2005 is compared to a model spectrum where 
only C II and O II line opacity has been included. The top panel is for solar abundances 
and the bottom panel using CNO enhanced abundances. The observed spectrum has been 
dereddened using E(B — V) — 0.035. 
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5. Day 17 

To examine the evolution of SN 2005cs we studied the spectrum obtained on July 14, 
2005. Using our results obtained above we will take the fiducial reddening to be E(B — V) — 
0.035. Figure [9] (top panel) shows our best fit model with solar abundances. The parameters 
for this model were: T mode i = 6000 K, v = 4000 km s _1 , and n — 8. The overall fit is 
quite good. The only noticeable mismatch is the trough just bluewa rd of 5700 A is not wel l 



fit. The red part of this trough is almost certainly due to Na I D. iPastorello et al.l (120061 ) 
identify the feature just to the blue as due to Ba II A5854, which won't be present in our 
models without enhanced s-process abundances or considering NLTE effects of Ba II, which 
is beyond the scope of this work (see below). The bottom panel of Figure M shows the fit 
with the same model parameters but CNO enhanced abundances. There is no noticeable 
degradation in the fit, and the enhanced sodium from CNO processing leads to a feature in 
the right place to fit the red part of the feature near 6000 A. To further test the sensitivity 
to the sodium abundance, we calculated a model with solar abundances, but with sodium 
enhanced by a factor of ten by number. Keeping the other parameters fixed, this model is 
shown in Figure [101 Clearly the enhanced sodium leads to a feature at the right place that 
is in fact too strong. Thus, there is evidence for enhanced sodium which would most likely 
come from CNO processing. Overall the sodium enhanced model is a bit too cool, but we 
have held the other parameters fixed in order to illustrate just the effect of enhancing the 
sodium abundance. 
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6. Day 34 

As our final epoch, we model spectra obtained on July 31, 2005. The parameters for 
this model were: T mo d e i = 5500 K, t> = 6000 km s _1 , and n = 12. This parameter set is 
somewhat odd, considering that it implies that the photosphere has actually moved outward 
from 4000 km s _1 on day 17 to 6000 km s _1 on day 31 and that the density profile is different. 
This contradiction indicates the limits of our modeling with a simple single power-law and 
uniform compositions. Figure [11] (top panel) shows the fit with solar abundances. Overall 
the color is well reproduced, Ha, and both the Ca II H+K and IR triplet features are well fit, 
but many other lines are not, including other lines in the Balmer series. In fact while the Ha 
absorption trough is quite well fit, the emission peak is too wide and too red. This is almost 
certainly due to our simple single power-law density profile as we will discuss in future work. 
Figure [11] (bottom panel) shows the fit with CNO processed abundances. As was the case 
for day 17, the only real improvement is that now Na I D appears in the synthetic spectrum, 
although it is too weak. Figure [T2] shows the result of enhancing Na by a factor of ten in 
the solar abundance mixture, now Na I D shows up nicely, although not quite at the right 
velocity or strength. 

Figure [TBI shows a model with solar abundances, T mo d e i = 5500 K, t> = 3000 km s _1 , 
and n = 12. While narrow lines appear, neither the Balmer lines nor the calcium lines are 
well fit and none of the absorption troughs appear at the right velocity. 



Pastorello et all (120061 ) found evidence for s-process elements such as barium, and charged 
particle elements such as scandium, which if present in our models would certainly improve 
the fit in Fig. [TT] In our models these species are treated in LTE, but it is possible that 
enhanced abundances would be required to produce th e strong lines observed. In studying 



the SN 1987A analog SN 1998A iPastorello et all (120051 ) concluded that the barium lines are 



quite temperature sensitive and thus the difference in strength between these lines and the 



very s trong lines observed in 1998A could well be due to NLTE effects. lUtrobin fc Chugai 



( 120051 ) also found that the barium lines in SN 1987A was strongly dependent on the radi- 
ation field and that a full NLTE solution with a reduced UV flux due to line-blanketing 
reduces the need for enhanced barium abundance. Thus, with a better model for this epoch 
and treating barium and scandium in NLTE we are likely to get significantly improved fits, 
without required enhanced s-process abundances. This is a subject for future work. 
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Fig. 9. — . The observed spectrum on Jul 14, 2005 is compared to a model spectrum using 
solar abundances (top panel) and CNO enhanced abundances (bottom panel). The observed 
spectrum has been dereddened using E(B — V) = 0.035. 
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Fig. 10. — . The observed spectrum on Jul 14, 2005 is compared to a model spectrum using 
solar abundances but sodium has been enhanced by a factor of 10. The observed spectrum 
has been dereddened using E(B — V) = 0.035. 
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Fig. 11. — . The observed spectrum on Jul 31, 2005 is compared to a model spectrum using 
solar abundances (top panel) and CNO enhanced abundances (bottom panel). The observed 
spectrum has been dereddened using E(B — V) = 0.035. 
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Fig. 12. — . The observed spectrum on Jul 31, 2005 is compared to a model spectrum using 
solar abundances but sodium has been enhanced by a factor of 10. The observed spectrum 
has been dereddened using E(B — V) = 0.035. 
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Fig. 13. — . The observed spectrum on Jul 31, 2005 is compared to a model spectrum using 
solar abundances and v = 3000 km s" 1 . The observed spectrum has been dereddened using 
E(B -V) = 0.035. 
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Distance 



One method of deter mining distances using 



ype II supernovae i s the "expanding pho- 



tosphere method" 


5PM, 


Kirshner & Kwan 


1974; 


Branch et a 


. 1981 


; Eastman & Kirshner 


1989; 


Schmidt et al. 


1994; 


Eastman et al. 


1996: Irlamuv et al. 


2001; : 


Leonard et al.l 


2002) a 


variation of the Baade-Wesselink method ( 


Baade Il926 


). The EPM method assumes that for 



SNe IIP, with intact hydrogen envelopes, the spectrum is not far from that of a blackbody 
and hence the luminosity is approximately given by 

L = 4vr C 2 R 2 o T 4 

where R is the radius of the photosphere, T is the effective temperature, a is the radia- 
tion constant, and ( is the "dilution factor" whi ch takes into account that in a scatteri ng 
dominated atmosphere the bla ckbody is diluted (jHershkowitz. Linder. &; Wagonerlll986al Jb: 
Hershkowitz fc Wagonerl 119871 ). The temperature is found from observed colors, so in fact 
is a color temperature and not an effective temperature, the photospheric velocity can be 
estimated from observed spectra using the velocities of the weakest lines, 



R = vt, 



the dilution factor is estimated from synthetic spectral models, and t comes from the light 
curve and demanding self-consistency. 

Both an advantage and disadvantage of EPM is that it primarily requires photometry. 
Spectra are only used to determine the photospheric velocity, colors yield the color tem- 
perature, which in turn is used to determine the appropriate dilution factor (from model 
results). This method suffers from uncertainties in determining the dilution factors, the 
difficulty of knowing which lines to use as velocity indicators, uncertainties between color 
temperatures and effective temperatures, and questions of how to match the photospheric 
radius used in the models to determine the dilution factor and the radius of the line form- 



ing region (IHamuy et al.l 1200 It iLeonard et al.l 120021) . In spite of this the EPM method was 
successfully applied to SN 1987A in the LMC faastman fc Kirshnerl Il989l : branch! Il987h 
which led to hopes that the EPM method would lead to accurate distances, independent 
of other astronomical cali brators. Recently, the EPM method was applied to the very we ll 



observed SN IIP 1999em ( Hamuy et al.ll200ll ; ILeonard et al. 



2002 



Elmhamdi et al.l 120031 ). 



All three groups found a distance of 7.5-8.0 Mpc. ILeonard et al.l (120031 ) subsequently used 
HST to obtain a Cepheid distance to the parent galaxy of SN 1999em, NGC 1637, and 
found 11.7 ± 1.0 Mpc, a value 50% larger than that obtained with EPM. iDessart fc Hillier 
( 120061 ) calculated updated values for the dilution factor and found good agreement with the 
Cepheid distance using EPM, but emphasized that dilution factors would have to be calcu- 
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lated for each supernova individually. Using SEAM iBaron et al.l (120041 ) found a distance to 
SN 1999em of 12.5 ± 2.3 Mpc. 

With modern detailed NLTE radiative transfer codes, accurate synthetic spectra of 
all types of supe rnovae can be calculated . The Spectral-fitting Expanding Atm osphere 
Method (SEAM, baron et al1ll995al . Il996l : iLentz et al.lboOll : iMitchell et aD l2002h was de- 
veloped using the generali zed stellar atmosphere code PHOENIX (for a review of the code see 
Hauschildt fc Baron! Il999i ). While SEAM is similar to EPM in spirit, it avoids the use of 
dilution factors and color temperatures. Velocities are determined accurately by actually 
fitting synthetic and observed spectra. The radius is still determined by the relationship 
R = vt, (which is an excellent approximation because all supernovae quickly reach homolo- 
gous expansion) and the explosion time is found by demanding self consistency. SEAM uses 
all the spectral information available in the observed spectra simultaneously which broad- 
ens the base of parameter determination. Since the spectral energy distribution is known 
completely from the calculated synthetic spectra, one may calculate the absolute magnitude, 
Mx, in any photometric band X, 



M 



x 



poo 

-2.5 log / Sx(\)L x d\ + C x 
Jo 



where Sx is the response of filter X, L\ is the luminosity per unit wavelength, and Cx is 
the zero point of filter X determined from standard stars. Then one immediately obtains a 
distance modulus fix, which is a measure of the distance 

A*x = m x - M x - A x - K x = 5 log (d/lOpc), 

where m x is the apparent magnitude in band X, A x is the extinction due to dust along 
the line of sight both in the host galaxy and in our own galaxy, and Kx is the k-correction 
( Oke fc Sandagel 1968 ). The SEAM method does not need to invoke a blackbody assumption 
or to calculate dilution factors. 

We calculate a SEAM distance to SN 2005cs using our first two early dates, day 5 and 
day 17, since our best fit to day 34 was contradictory in that the photosphere moved out 
and the density slope was different at a previous time. With only two day s we don't have 



much leverage on the time of explosion, thus we adopt the explosion time of iPastorello et al. 



( 120061 ). Using 5 photometric bands UBVRI and our adopted reddening of E(B — V) = 0.035 
we find /i = 29.5 ± 0.18 which corresponds to a distance of 7.9lgg Mpc where the error is 
just the formal 1 — a e rror and does not include systematic error. Using the EPM method 
akats &: Vinkol (120061) found D = 7.59 ± 1.02 Mpc, and using a standard candle method 



( IHamuy 



2003 



Nugent et al.ll2006f ) they found D = 6.02±1.3 Mpc. Using the S EAM method 



on SN 19941 a Type Ic that is also in M51 we found D = 6.02 ± 1.92 Mpc (IBaron et al. 
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19961 ). Given that SNe Ic are much more difficult to model than SNe II, the fact that the two 



SEAM results agree within the errors is strong evidence that not only is the SE AM method 
robust, but the errors can be reasonably well estimated. iTakats fe Vinkol (120061 ) summarize 
the other distance estimates obtained to M51 and find an average value of D — 7.1±1.2 Mpc. 
Given the wide range of distances obtained and the relative closeness of M51, an effort should 
be made to obtain a Cepheid distance to M51. 



8. Discussion and Conclusions 

We have calculated a small grid of simple parameterized explosion models for SN 2005cs. 
We have shown that the reddening estimated in previous studies is likely too high and we 
adopted the lowest possible value, that of foreground reddening alone. Additionally we 
have shown that while there is evidence for enhanced nitrogen and sodium as would be 
expected from CNO processing in the pre-supernova envelope, the specific set of abundances 
that we adopted were too depleted in oxygen to match the observed spectrum. We have not 
attempted to exhaust the full parameter space and determine exact abundances in this work, 
merely to show that such an approach is feasible given sufficient computational resources. 
Additionally, we had trouble obtaining a physically plausible fit to the final epoch we studied, 
34 days past explosion. All of these shortcomings point that a better approach would be to 
run large grids using the best current stellar evolution models and compare the results to 
many well-observed SNe IIP. This should be feasible in the near future. Using the SEAM 
method we obtained a distance to M51 which is somewhat larger than the accepted distance. 
The SEAM method has been shown to be reliable in modeling SNe IIP and our result is 
reliable enough that it should be calibrated with a Cepheid distance to M51. This could 
then be compared with the results from a larger more homogeneous grid. 
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